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The methods of light-front quantization and Pauli-Villars regularization are applied to a nonper- 
turbative calculation of the dressed-electron state in quantum electrodynamics. This is intended 
as a test of the methods in a gauge theory, as a precursor to possible methods for the nonperturba- 
tive solution of quantum chromodynamics. The electron state is truncated to include at most two 
photons and no positrons in the Fock basis, and the wave functions of the dressed state are used to 
compute the electrons's anomalous magnetic moment. A choice of regularization that preserves 
the chiral symmetry of the massless limit is critical for the success of the calculation. 
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1. Introduction 

The purpose of this work is to explore a nonperturbative method that can be used to solve for 
the bound states of quantum field theories, in particular QCD. The problem is notoriously diffi- 
cult, and there are only a few approaches. These include lattice gauge theory [jlj], the transverse 
lattice 0], Dyson-Schwinger equations [||], Bethe-Salpeter equation, similarity transformations 
combined with construction of effective fields JjJ], and light-front Hamiltonians with either stan- 
dard [Q] or sector-dependent parameterizations [0, [7L [8p. We use the light-front Hamiltonian ap- 
proach with Pauli-Villars (PV) [^] regularizaton and standard parameterization, where the bare 
parameters of the Lagrangian do not depend on the Fock sector. As a test in a gauge theory, we 
consider light-front QED and specifically the eigenstate of the dressed electron and its anomalous 
moment (fij US 01- 



We use light-cone coordinates [14, 15], chosen in order to have well-defined Fock-state ex- 
pansions and a simple vacuum. The time coordinate is x + =t + z and the space coordinates are 
x = (x~~ ,x±), with x~ = t — z and x± = (x,y). The light-cone energy is p" = E — p z , and the three- 
momentum is p = (p + ,p±), with p + = E + p z and p± = (p x ,p y ). The mass-shell condition p 2 = m 2 

becomes p~~ = - p + ± . The simple vacuum follows from the positivity of the plus component of 



the momentum: p + = ym 2 + p 2 + p\+ p z > 0. 

To regulate QED, we use the Pauli-Villars technique The basic idea is to subtract from 
each integral a contribution of the same form but of a PV particle with a much larger mass. This 
can be done by adding negative metric particles to the Lagrangian. A particular advantage of 
PV regularization is preservation of at least some symmetries; in particular, it is automatically 
relativistically covariant. 

From the PV-regulated light-front QED Lagrangian, we construct the Hamiltonian and 
solve the mass eigenvalue problem <^~|P) = j+\P) in the approximation that the electron eigen- 
state is a truncated Fock-state expansion with at most two photons and no positrons. From this 
approximate eigenstate, we compute the anomalous magnetic moment, as a test of the method. 

2. Light-front QED 

The light-front QED Lagrangian with one PV fermion and two PV photons is 

i=o L ^ z z 

i 

+ L ( ~ 1 )' W (»V - mi) Yi ~ e xj/'f yA^ , 

(=0 



(2.1) 



with 



1 2 

\jf = £ y/i, A M = £ \/%A ¥ , F iflv = d^Aiv - d v A^ . (2.2) 



i=0 i=0 

-2 



The coupling coefficients ^ are constrained by £,q = 1 and L/=o( — l)'^' = 0, and the requirement 
of chiral symmetry restoration in the limit of zero electron mass. At one loop, the chiral symmetry 
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2/2 

constraint becomes [1C] Ya= o(~ ^ ™i Y % m (Mf / m i) = 0; for nonperturbative solutions with 



more than one photon in the basis, the constraint must be imposed numerically [ 13] 



The light-front Hamiltonian without antifermion terms is then of the form [12] 



2 2 

^ -LI dA^(-i)%{p)bUp) 



+1 

l,H 



+ I 



(-l)'e%(*Hx(*) 

dpdq I fet, (>) b jjS (q) V? 2s (p, q) 



(2.3) 



with £^ = (—1,1,1,1) and the vertex functions given in JT^]. 

We work in a frame where the total transverse momentum P± is zero and expand the eigen- 
f unction for the dressed-fermion state with total J 7 = ±A in a Fock basis as 



w ± (D) = j>Lcaio> + 1 / dkc£(k)bi(p-k)"U® i°) 



(2.4) 



where we have truncated the expansion to include at most two photons. The z,- are the amplitudes 
for the bare electron states, with i = for the physical electron and i = 1 for the PV electron. 
The Cfjf are the two-body wave functions for Fock states with an electron of flavor i and spin 
component s and a photon of flavor j = 0, 1 or 2 and field component pi, expressed as functions 
of the photon momentum. The upper index of ± refers to the J z value of ±^ for the eigenstate. 
Similarly, the C-J fe are the three-body wave functions for the states with one electron and two 
photons, with flavors j and k and field components p: and v. 

The Fock expansion is an eigenstate of the light-front Hamiltonian with eigenvalue 
M 2 /P + . The wave functions then satisfy the following coupled integral equations: 



[M 2 -mj] Zi 



dq £ 

j,l,H 



(2.5) 



U?UP-q,P)c£(q) 



M 



2 m ]+q\ tf+q\ 



(i-y) y 



(2.6) 



2^ 



+ £(_ 1 )«+* e v d ! 



V:*{P-i-q,P-i)Cl£U,q) 
+U2 (P-q~q,P- q)CZt (f, g) 
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M 



2 mj + (q l± + q 2± ) 2 Hj+q 1± ^lf+q 2 2± 



[l -yi -y 2 ) 



>'2 



Q 



(2.7) 



+U?_ s {P-q l -q 2 ,P- q 2 )C al _ s ^ 
+ Vb[Vte(P-q 1 -q 2 ,P-q 1 )C%( Sl ) 

+ Uy a! . s (P-q, -g 2 ,P-g 1 )C^(£ 1 )] } . 

The anomalous moment a e can be computed from the spin-flip matrix element of the electro 
magnetic current /+ [|T6|]. At zero momentum transfer, we have a e = F 2 (0) and 

/ i ;/2+l /_ lV "+i' \* 



(2.8) 



7=0,2 



xy[ dk x +l d£ 



K i'=0k'=j/2 V W 
K i=0k=j/2 VW , 



sfiv' j,k=0,2 

x ft j L k L ( - 1) '" +r+m ' ^ c ^ ( ^^ ), 

\i'=QV=j/2m'=k/2 VZl'Zm' a/1 + 5/' m '5 MV 
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'tie *iT (- i )'' +/+m ^ c "lr(^^) 

v i=0 l=j/2 m=k/2 Vwlm a/ 1 + Wfiv 



The terms that depend on the three-body wave functions C^f are higher order in a than the leading 
two-body terms. Given the numerical errors in the leading terms, these three-body contributions 
are not significant and are not evaluated. The important three-body contributions come from the 
couplings of the three-body wave functions that will enter the calculation of the two-body wave 
functions. 



3. Solution of the Coupled Equations 



The first and third equations of the coupled system, ( |2.5[ ) and ( |2.7| ), can be solved for the 
bare-electron amplitudes and one-electron/two-photon wave functions, respectively, in terms of the 
one-electron/one-photon wave functions. Substitution of these solutions into the second integral 



equation (2.6) yields a reduced integral eigenvalue problem in the one-electron/one-photon sector: 



M 



2 mj + q 2 ± Hj+qi 
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There is a total of 48 coupled equations, with i = 0,1; j = 0, 1,2; s = rb^; and ;U = ±, (±). 
The first term on the right-hand side of (3.1) is the self-energy contribution Jl2|]: 



1 dx d 2 k ± m i m i' ~ 2 -[=r m a + ■ 



a,o /v/ i — jc * 

with 

A/ ^ 2 + (1-^-^1 (3.3) 
The kernels and 7^ 2 ^ in the second and third terms correspond to interactions with zero or two 



photons in intermediate states. Details of these kernels can be found in [ 13] and [11]. 



The presence of the flavor changing self-energies, the /,-/,■/ with i ^ i', generates a fermion flavor 



mixing of the two-body wave functions []12|]. To resolve this, we write the integral equations for 
these wave functions in the form 

where A, ; - and Bj are defined by 

A d±A + tilA + ^JiJL_ M 2 i B c L hj^ = _o L Joj L 
1 l-y y 2%\-y 3 2%\-y 2k l-y 

and jf- is given by 



€ = I £V / w^£ ) ^v(^^;/,^)c;jr s ,(y,^) (3.6) 

i'j's'v Jo 

£ £ V I dy'dqfjjj^J^i^,q ± ;y',q 1 ^)Cy,f s ,(j',q_ L ). 

,■/;/„/„ JO 



i'j's'v 

We then construct wave functions that mix fermion flavors and diagonalize the left-hand side of 
Q: ygf = A tj cff s + (- 1 YBjC^ js . Solution of the resulting integral equations for the jfgf 
yields a as a function of mo and the PV masses. Then for given values of PV masses, we can seek 
the value of mo for which a takes the standard physical value e 2 /4n. The eigenproblem solution 
also yields the functions ffj s which determine the wave functions C^f. From these wave functions 
we can compute physical quantities as expectation values with respect to the projection [|l3|] of the 
eigenstate onto the physical subspace. 

The eigenvalue problem must first be solved for M = 0, with the coupling strength parameter 
(§2 adjusted to yield mo = 0. This determines the value of ^ that restores the chiral limit nonpertur- 
batively. The eigenvalue problem can then be solved for M = m e , the physical mass of the electron, 
and the anomalous moment calculated. 

If we retain only the self-energy contributions from the two-photon intermediate states, the 
equations for the two-body wave functions become much simpler, and the coupled integral equa- 
tions can be reduced to the one-electron sector. There, they can be solved analytically, except for 



the calculation of certain integrals [12]. 
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Figure 1: The anomalous moment of the electron in units of the Schwinger term (a/2n) plotted versus the 
PV photon mass, jUi, with the second PV photon mass, ji2, set to \/2p.\ and the PV electron mass m\ equal to 
2 • 10 4 »v The solid squares are the result of the full two-photon truncation with the correct, nonperturbative 
chiral constraint [Q. The open squares come from use of a perturbative, one-loop constraint. Results for 
the one-photon truncation [ |Io| ] (solid line) and the one-photon truncation with the two-photon self-energy 
contribution Jl2] ] (filled circles) are included for comparison. The resolutions used for the two-photon results 
are K — 50 to 150, combined with extrapolation to K = °o, and N± = 20. 



4. Results 

From the solutions to the eigenvalue problems, we compute the anomalous moment at fixed PV 
masses and fixed numerical resolution. We then study the behavior first as a function of the numer- 
ical resolution, which requires extrapolation, and then as a function of PV masses. The numerical 
resolution is marked by two parameters, K and N±, which control the number of quadrature points 
used in the longitudinal and transverse directions. The numerical convergence and extrapolation 
are illustrated in [O]. 

The results of the extrapolations are plotted in Fig. |]. Each value is close to the standard 
Schwinger result of a/2n and independent of }X\, to within numerical error. The results with only 
the two-photon self-energy contribution are actually better than the full two-photon results. This 
discrepancy should be due to the absence of electron-positron contributions, which are of the same 
order in a as the two-photon contributions; without the electron-positron contributions, we lack 
the cancellations that typically take place between contributions of the same order. 

We also see that the inclusion of the self-energy contribution is a significant improvement over 
the one-photon truncation. Thus, we expect that inclusion of three-photon self-energy effects will 
improve the two-photon results. 

Figure [j] also includes results obtained for the two-photon truncation when only the one-loop 
chiral constraint is satisfied. Without the full nonperturbative constraint, the results are very sensi- 
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tive to the PV photon mass pL\. This behavior repeats the pattern observed in p0| ] for a one-photon 
truncation without the corresponding one-loop constraint. The resulting /ii dependence is illus- 



trated in Fig. 2 of [JlOp. Thus, a successful calculation requires that the symmetry of the chiral limit 
be maintained. 
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